


<!DOCTYPE html>
<!--[if IE 8]><html class="no-js lt-ie9" lang="en" > <![endif]-->
<!--[if gt IE 8]><!--> <html class="no-js" lang="en" > <!--<![endif]-->
<head>
  <meta charset="utf-8">
  
  <meta name="viewport" content="width=device-width, initial-scale=1.0">
  
  <title>torch.distributions.kl &mdash; PyTorch master documentation</title>
  

  
  
  
  
    <link rel="canonical" href="https://pytorch.org/docs/stable/_modules/torch/distributions/kl.html"/>
  

  

  
  
    

  

  <link rel="stylesheet" href="../../../_static/css/theme.css" type="text/css" />
  <!-- <link rel="stylesheet" href="../../../_static/pygments.css" type="text/css" /> -->
  <link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/katex@0.10.0-beta/dist/katex.min.css" type="text/css" />
  <link rel="stylesheet" href="../../../_static/css/jit.css" type="text/css" />
  <link rel="stylesheet" href="https://cdn.jsdelivr.net/npm/katex@0.11.1/dist/katex.min.css" type="text/css" />
  <link rel="stylesheet" href="../../../_static/katex-math.css" type="text/css" />
    <link rel="index" title="Index" href="../../../genindex.html" />
    <link rel="search" title="Search" href="../../../search.html" /> 

  
  <script src="../../../_static/js/modernizr.min.js"></script>

  <!-- Preload the theme fonts -->

<link rel="preload" href="../../../_static/fonts/FreightSans/freight-sans-book.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="../../../_static/fonts/FreightSans/freight-sans-medium.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="../../../_static/fonts/IBMPlexMono/IBMPlexMono-Medium.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="../../../_static/fonts/FreightSans/freight-sans-bold.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="../../../_static/fonts/FreightSans/freight-sans-medium-italic.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="../../../_static/fonts/IBMPlexMono/IBMPlexMono-SemiBold.woff2" as="font" type="font/woff2" crossorigin="anonymous">

<!-- Preload the katex fonts -->

<link rel="preload" href="https://cdn.jsdelivr.net/npm/katex@0.10.0/dist/fonts/KaTeX_Math-Italic.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="https://cdn.jsdelivr.net/npm/katex@0.10.0/dist/fonts/KaTeX_Main-Regular.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="https://cdn.jsdelivr.net/npm/katex@0.10.0/dist/fonts/KaTeX_Main-Bold.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="https://cdn.jsdelivr.net/npm/katex@0.10.0/dist/fonts/KaTeX_Size1-Regular.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="https://cdn.jsdelivr.net/npm/katex@0.10.0/dist/fonts/KaTeX_Size4-Regular.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="https://cdn.jsdelivr.net/npm/katex@0.10.0/dist/fonts/KaTeX_Size2-Regular.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="https://cdn.jsdelivr.net/npm/katex@0.10.0/dist/fonts/KaTeX_Size3-Regular.woff2" as="font" type="font/woff2" crossorigin="anonymous">
<link rel="preload" href="https://cdn.jsdelivr.net/npm/katex@0.10.0/dist/fonts/KaTeX_Caligraphic-Regular.woff2" as="font" type="font/woff2" crossorigin="anonymous">
</head>

<div class="container-fluid header-holder tutorials-header" id="header-holder">
  <div class="container">
    <div class="header-container">
      <a class="header-logo" href="https://pytorch.org/" aria-label="PyTorch"></a>

      <div class="main-menu">
        <ul>
          <li>
            <a href="https://pytorch.org/get-started">Get Started</a>
          </li>

          <li>
            <div class="ecosystem-dropdown">
              <a id="dropdownMenuButton" data-toggle="ecosystem-dropdown">
                Ecosystem
              </a>
              <div class="ecosystem-dropdown-menu">
                <a class="nav-dropdown-item" href="https://pytorch.org/hub"">
                  <span class=dropdown-title>Models (Beta)</span>
                  <p>Discover, publish, and reuse pre-trained models</p>
                </a>
                <a class="nav-dropdown-item" href="https://pytorch.org/ecosystem">
                  <span class=dropdown-title>Tools & Libraries</span>
                  <p>Explore the ecosystem of tools and libraries</p>
                </a>
              </div>
            </div>
          </li>

          <li>
            <a href="https://pytorch.org/mobile">Mobile</a>
          </li>

          <li>
            <a href="https://pytorch.org/blog/">Blog</a>
          </li>

          <li>
            <a href="https://pytorch.org/tutorials">Tutorials</a>
          </li>

          <li class="active">
            <a href="https://pytorch.org/docs/stable/index.html">Docs</a>
          </li>

          <li>
            <div class="resources-dropdown">
              <a id="resourcesDropdownButton" data-toggle="resources-dropdown">
                Resources
              </a>
              <div class="resources-dropdown-menu">
                <a class="nav-dropdown-item" href="https://pytorch.org/resources"">
                  <span class=dropdown-title>Developer Resources</span>
                  <p>Find resources and get questions answered</p>
                </a>
                <a class="nav-dropdown-item" href="https://pytorch.org/features">
                  <span class=dropdown-title>About</span>
                  <p>Learn about PyTorch’s features and capabilities</p>
                </a>
              </div>
            </div>
          </li>

          <li>
            <a href="https://github.com/pytorch/pytorch">Github</a>
          </li>
        </ul>
      </div>

      <a class="main-menu-open-button" href="#" data-behavior="open-mobile-menu"></a>
    </div>

  </div>
</div>


<body class="pytorch-body">

   

    

    <div class="table-of-contents-link-wrapper">
      <span>Table of Contents</span>
      <a href="#" class="toggle-table-of-contents" data-behavior="toggle-table-of-contents"></a>
    </div>

    <nav data-toggle="wy-nav-shift" class="pytorch-left-menu" id="pytorch-left-menu">
      <div class="pytorch-side-scroll">
        <div class="pytorch-menu pytorch-menu-vertical" data-spy="affix" role="navigation" aria-label="main navigation">
          <div class="pytorch-left-menu-search">
            

            
              
              
                <div class="version">
                  master (1.5.0 )
                </div>
              
            

            


  


<div role="search">
  <form id="rtd-search-form" class="wy-form" action="../../../search.html" method="get">
    <input type="text" name="q" placeholder="Search Docs" />
    <input type="hidden" name="check_keywords" value="yes" />
    <input type="hidden" name="area" value="default" />
  </form>
</div>

            
          </div>

          
<div>
  <a style="color:#F05732" href="https://pytorch.org/docs/stable/_modules/torch/distributions/kl.html">
    You are viewing unstable developer preview docs.
    Click here to view docs for latest stable release.
  </a>
</div>

            
            
              
            
            
              <p class="caption"><span class="caption-text">Notes</span></p>
<ul>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/amp_examples.html">Automatic Mixed Precision examples</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/autograd.html">Autograd mechanics</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/broadcasting.html">Broadcasting semantics</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/cpu_threading_torchscript_inference.html">CPU threading and TorchScript inference</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/cuda.html">CUDA semantics</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/ddp.html">Distributed Data Parallel</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/extending.html">Extending PyTorch</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/faq.html">Frequently Asked Questions</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/large_scale_deployments.html">Features for large-scale deployments</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/multiprocessing.html">Multiprocessing best practices</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/randomness.html">Reproducibility</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/serialization.html">Serialization semantics</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../notes/windows.html">Windows FAQ</a></li>
</ul>
<p class="caption"><span class="caption-text">Language Bindings</span></p>
<ul>
<li class="toctree-l1"><a class="reference external" href="https://pytorch.org/cppdocs/">C++ API</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../packages.html">Javadoc</a></li>
</ul>
<p class="caption"><span class="caption-text">Python API</span></p>
<ul>
<li class="toctree-l1"><a class="reference internal" href="../../../torch.html">torch</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../nn.html">torch.nn</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../nn.functional.html">torch.nn.functional</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../tensors.html">torch.Tensor</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../tensor_attributes.html">Tensor Attributes</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../tensor_view.html">Tensor Views</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../autograd.html">torch.autograd</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../cuda.html">torch.cuda</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../amp.html">torch.cuda.amp</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../distributed.html">torch.distributed</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../distributions.html">torch.distributions</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../hub.html">torch.hub</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../jit.html">torch.jit</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../nn.init.html">torch.nn.init</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../onnx.html">torch.onnx</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../optim.html">torch.optim</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../quantization.html">Quantization</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../rpc/index.html">Distributed RPC Framework</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../random.html">torch.random</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../sparse.html">torch.sparse</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../storage.html">torch.Storage</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../bottleneck.html">torch.utils.bottleneck</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../checkpoint.html">torch.utils.checkpoint</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../cpp_extension.html">torch.utils.cpp_extension</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../data.html">torch.utils.data</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../dlpack.html">torch.utils.dlpack</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../model_zoo.html">torch.utils.model_zoo</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../tensorboard.html">torch.utils.tensorboard</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../type_info.html">Type Info</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../named_tensor.html">Named Tensors</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../name_inference.html">Named Tensors operator coverage</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../__config__.html">torch.__config__</a></li>
</ul>
<p class="caption"><span class="caption-text">Libraries</span></p>
<ul>
<li class="toctree-l1"><a class="reference external" href="https://pytorch.org/audio">torchaudio</a></li>
<li class="toctree-l1"><a class="reference external" href="https://pytorch.org/text">torchtext</a></li>
<li class="toctree-l1"><a class="reference external" href="https://pytorch.org/elastic/">TorchElastic</a></li>
<li class="toctree-l1"><a class="reference external" href="https://pytorch.org/serve">TorchServe</a></li>
<li class="toctree-l1"><a class="reference external" href="http://pytorch.org/xla/">PyTorch on XLA Devices</a></li>
</ul>
<p class="caption"><span class="caption-text">Community</span></p>
<ul>
<li class="toctree-l1"><a class="reference internal" href="../../../community/contribution_guide.html">PyTorch Contribution Guide</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../community/governance.html">PyTorch Governance</a></li>
<li class="toctree-l1"><a class="reference internal" href="../../../community/persons_of_interest.html">PyTorch Governance | Persons of Interest</a></li>
</ul>

            
          

        </div>
      </div>
    </nav>

    <div class="pytorch-container">
      <div class="pytorch-page-level-bar" id="pytorch-page-level-bar">
        <div class="pytorch-breadcrumbs-wrapper">
          















<div role="navigation" aria-label="breadcrumbs navigation">

  <ul class="pytorch-breadcrumbs">
    
      <li>
        <a href="../../../index.html">
          
            Docs
          
        </a> &gt;
      </li>

        
          <li><a href="../../index.html">Module code</a> &gt;</li>
        
          <li><a href="../../torch.html">torch</a> &gt;</li>
        
      <li>torch.distributions.kl</li>
    
    
      <li class="pytorch-breadcrumbs-aside">
        
      </li>
    
  </ul>

  
</div>
        </div>

        <div class="pytorch-shortcuts-wrapper" id="pytorch-shortcuts-wrapper">
          Shortcuts
        </div>
      </div>

      <section data-toggle="wy-nav-shift" id="pytorch-content-wrap" class="pytorch-content-wrap">
        <div class="pytorch-content-left">

        
          
          <div class="rst-content">
          
            <div role="main" class="main-content" itemscope="itemscope" itemtype="http://schema.org/Article">
             <article itemprop="articleBody" id="pytorch-article" class="pytorch-article">
              
  <h1>Source code for torch.distributions.kl</h1><div class="highlight"><pre>
<span></span><span class="kn">import</span> <span class="nn">math</span>
<span class="kn">import</span> <span class="nn">warnings</span>
<span class="kn">from</span> <span class="nn">functools</span> <span class="kn">import</span> <span class="n">total_ordering</span>

<span class="kn">import</span> <span class="nn">torch</span>
<span class="kn">from</span> <span class="nn">torch._six</span> <span class="kn">import</span> <span class="n">inf</span>

<span class="kn">from</span> <span class="nn">.bernoulli</span> <span class="kn">import</span> <span class="n">Bernoulli</span>
<span class="kn">from</span> <span class="nn">.beta</span> <span class="kn">import</span> <span class="n">Beta</span>
<span class="kn">from</span> <span class="nn">.binomial</span> <span class="kn">import</span> <span class="n">Binomial</span>
<span class="kn">from</span> <span class="nn">.categorical</span> <span class="kn">import</span> <span class="n">Categorical</span>
<span class="kn">from</span> <span class="nn">.continuous_bernoulli</span> <span class="kn">import</span> <span class="n">ContinuousBernoulli</span>
<span class="kn">from</span> <span class="nn">.dirichlet</span> <span class="kn">import</span> <span class="n">Dirichlet</span>
<span class="kn">from</span> <span class="nn">.distribution</span> <span class="kn">import</span> <span class="n">Distribution</span>
<span class="kn">from</span> <span class="nn">.exponential</span> <span class="kn">import</span> <span class="n">Exponential</span>
<span class="kn">from</span> <span class="nn">.exp_family</span> <span class="kn">import</span> <span class="n">ExponentialFamily</span>
<span class="kn">from</span> <span class="nn">.gamma</span> <span class="kn">import</span> <span class="n">Gamma</span>
<span class="kn">from</span> <span class="nn">.geometric</span> <span class="kn">import</span> <span class="n">Geometric</span>
<span class="kn">from</span> <span class="nn">.gumbel</span> <span class="kn">import</span> <span class="n">Gumbel</span>
<span class="kn">from</span> <span class="nn">.half_normal</span> <span class="kn">import</span> <span class="n">HalfNormal</span>
<span class="kn">from</span> <span class="nn">.independent</span> <span class="kn">import</span> <span class="n">Independent</span>
<span class="kn">from</span> <span class="nn">.laplace</span> <span class="kn">import</span> <span class="n">Laplace</span>
<span class="kn">from</span> <span class="nn">.lowrank_multivariate_normal</span> <span class="kn">import</span> <span class="p">(</span><span class="n">LowRankMultivariateNormal</span><span class="p">,</span> <span class="n">_batch_lowrank_logdet</span><span class="p">,</span>
                                          <span class="n">_batch_lowrank_mahalanobis</span><span class="p">)</span>
<span class="kn">from</span> <span class="nn">.multivariate_normal</span> <span class="kn">import</span> <span class="p">(</span><span class="n">MultivariateNormal</span><span class="p">,</span> <span class="n">_batch_mahalanobis</span><span class="p">)</span>
<span class="kn">from</span> <span class="nn">.normal</span> <span class="kn">import</span> <span class="n">Normal</span>
<span class="kn">from</span> <span class="nn">.one_hot_categorical</span> <span class="kn">import</span> <span class="n">OneHotCategorical</span>
<span class="kn">from</span> <span class="nn">.pareto</span> <span class="kn">import</span> <span class="n">Pareto</span>
<span class="kn">from</span> <span class="nn">.poisson</span> <span class="kn">import</span> <span class="n">Poisson</span>
<span class="kn">from</span> <span class="nn">.transformed_distribution</span> <span class="kn">import</span> <span class="n">TransformedDistribution</span>
<span class="kn">from</span> <span class="nn">.uniform</span> <span class="kn">import</span> <span class="n">Uniform</span>
<span class="kn">from</span> <span class="nn">.utils</span> <span class="kn">import</span> <span class="n">_sum_rightmost</span>

<span class="n">_KL_REGISTRY</span> <span class="o">=</span> <span class="p">{}</span>  <span class="c1"># Source of truth mapping a few general (type, type) pairs to functions.</span>
<span class="n">_KL_MEMOIZE</span> <span class="o">=</span> <span class="p">{}</span>  <span class="c1"># Memoized version mapping many specific (type, type) pairs to functions.</span>


<div class="viewcode-block" id="register_kl"><a class="viewcode-back" href="../../../distributions.html#torch.distributions.kl.register_kl">[docs]</a><span class="k">def</span> <span class="nf">register_kl</span><span class="p">(</span><span class="n">type_p</span><span class="p">,</span> <span class="n">type_q</span><span class="p">):</span>
    <span class="sd">&quot;&quot;&quot;</span>
<span class="sd">    Decorator to register a pairwise function with :meth:`kl_divergence`.</span>
<span class="sd">    Usage::</span>

<span class="sd">        @register_kl(Normal, Normal)</span>
<span class="sd">        def kl_normal_normal(p, q):</span>
<span class="sd">            # insert implementation here</span>

<span class="sd">    Lookup returns the most specific (type,type) match ordered by subclass. If</span>
<span class="sd">    the match is ambiguous, a `RuntimeWarning` is raised. For example to</span>
<span class="sd">    resolve the ambiguous situation::</span>

<span class="sd">        @register_kl(BaseP, DerivedQ)</span>
<span class="sd">        def kl_version1(p, q): ...</span>
<span class="sd">        @register_kl(DerivedP, BaseQ)</span>
<span class="sd">        def kl_version2(p, q): ...</span>

<span class="sd">    you should register a third most-specific implementation, e.g.::</span>

<span class="sd">        register_kl(DerivedP, DerivedQ)(kl_version1)  # Break the tie.</span>

<span class="sd">    Args:</span>
<span class="sd">        type_p (type): A subclass of :class:`~torch.distributions.Distribution`.</span>
<span class="sd">        type_q (type): A subclass of :class:`~torch.distributions.Distribution`.</span>
<span class="sd">    &quot;&quot;&quot;</span>
    <span class="k">if</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">type_p</span><span class="p">,</span> <span class="nb">type</span><span class="p">)</span> <span class="ow">and</span> <span class="nb">issubclass</span><span class="p">(</span><span class="n">type_p</span><span class="p">,</span> <span class="n">Distribution</span><span class="p">):</span>
        <span class="k">raise</span> <span class="ne">TypeError</span><span class="p">(</span><span class="s1">&#39;Expected type_p to be a Distribution subclass but got </span><span class="si">{}</span><span class="s1">&#39;</span><span class="o">.</span><span class="n">format</span><span class="p">(</span><span class="n">type_p</span><span class="p">))</span>
    <span class="k">if</span> <span class="ow">not</span> <span class="nb">isinstance</span><span class="p">(</span><span class="n">type_q</span><span class="p">,</span> <span class="nb">type</span><span class="p">)</span> <span class="ow">and</span> <span class="nb">issubclass</span><span class="p">(</span><span class="n">type_q</span><span class="p">,</span> <span class="n">Distribution</span><span class="p">):</span>
        <span class="k">raise</span> <span class="ne">TypeError</span><span class="p">(</span><span class="s1">&#39;Expected type_q to be a Distribution subclass but got </span><span class="si">{}</span><span class="s1">&#39;</span><span class="o">.</span><span class="n">format</span><span class="p">(</span><span class="n">type_q</span><span class="p">))</span>

    <span class="k">def</span> <span class="nf">decorator</span><span class="p">(</span><span class="n">fun</span><span class="p">):</span>
        <span class="n">_KL_REGISTRY</span><span class="p">[</span><span class="n">type_p</span><span class="p">,</span> <span class="n">type_q</span><span class="p">]</span> <span class="o">=</span> <span class="n">fun</span>
        <span class="n">_KL_MEMOIZE</span><span class="o">.</span><span class="n">clear</span><span class="p">()</span>  <span class="c1"># reset since lookup order may have changed</span>
        <span class="k">return</span> <span class="n">fun</span>

    <span class="k">return</span> <span class="n">decorator</span></div>


<span class="nd">@total_ordering</span>
<span class="k">class</span> <span class="nc">_Match</span><span class="p">(</span><span class="nb">object</span><span class="p">):</span>
    <span class="vm">__slots__</span> <span class="o">=</span> <span class="p">[</span><span class="s1">&#39;types&#39;</span><span class="p">]</span>

    <span class="k">def</span> <span class="fm">__init__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="o">*</span><span class="n">types</span><span class="p">):</span>
        <span class="bp">self</span><span class="o">.</span><span class="n">types</span> <span class="o">=</span> <span class="n">types</span>

    <span class="k">def</span> <span class="fm">__eq__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">other</span><span class="p">):</span>
        <span class="k">return</span> <span class="bp">self</span><span class="o">.</span><span class="n">types</span> <span class="o">==</span> <span class="n">other</span><span class="o">.</span><span class="n">types</span>

    <span class="k">def</span> <span class="fm">__le__</span><span class="p">(</span><span class="bp">self</span><span class="p">,</span> <span class="n">other</span><span class="p">):</span>
        <span class="k">for</span> <span class="n">x</span><span class="p">,</span> <span class="n">y</span> <span class="ow">in</span> <span class="nb">zip</span><span class="p">(</span><span class="bp">self</span><span class="o">.</span><span class="n">types</span><span class="p">,</span> <span class="n">other</span><span class="o">.</span><span class="n">types</span><span class="p">):</span>
            <span class="k">if</span> <span class="ow">not</span> <span class="nb">issubclass</span><span class="p">(</span><span class="n">x</span><span class="p">,</span> <span class="n">y</span><span class="p">):</span>
                <span class="k">return</span> <span class="kc">False</span>
            <span class="k">if</span> <span class="n">x</span> <span class="ow">is</span> <span class="ow">not</span> <span class="n">y</span><span class="p">:</span>
                <span class="k">break</span>
        <span class="k">return</span> <span class="kc">True</span>


<span class="k">def</span> <span class="nf">_dispatch_kl</span><span class="p">(</span><span class="n">type_p</span><span class="p">,</span> <span class="n">type_q</span><span class="p">):</span>
    <span class="sd">&quot;&quot;&quot;</span>
<span class="sd">    Find the most specific approximate match, assuming single inheritance.</span>
<span class="sd">    &quot;&quot;&quot;</span>
    <span class="n">matches</span> <span class="o">=</span> <span class="p">[(</span><span class="n">super_p</span><span class="p">,</span> <span class="n">super_q</span><span class="p">)</span> <span class="k">for</span> <span class="n">super_p</span><span class="p">,</span> <span class="n">super_q</span> <span class="ow">in</span> <span class="n">_KL_REGISTRY</span>
               <span class="k">if</span> <span class="nb">issubclass</span><span class="p">(</span><span class="n">type_p</span><span class="p">,</span> <span class="n">super_p</span><span class="p">)</span> <span class="ow">and</span> <span class="nb">issubclass</span><span class="p">(</span><span class="n">type_q</span><span class="p">,</span> <span class="n">super_q</span><span class="p">)]</span>
    <span class="k">if</span> <span class="ow">not</span> <span class="n">matches</span><span class="p">:</span>
        <span class="k">return</span> <span class="bp">NotImplemented</span>
    <span class="c1"># Check that the left- and right- lexicographic orders agree.</span>
    <span class="n">left_p</span><span class="p">,</span> <span class="n">left_q</span> <span class="o">=</span> <span class="nb">min</span><span class="p">(</span><span class="n">_Match</span><span class="p">(</span><span class="o">*</span><span class="n">m</span><span class="p">)</span> <span class="k">for</span> <span class="n">m</span> <span class="ow">in</span> <span class="n">matches</span><span class="p">)</span><span class="o">.</span><span class="n">types</span>
    <span class="n">right_q</span><span class="p">,</span> <span class="n">right_p</span> <span class="o">=</span> <span class="nb">min</span><span class="p">(</span><span class="n">_Match</span><span class="p">(</span><span class="o">*</span><span class="nb">reversed</span><span class="p">(</span><span class="n">m</span><span class="p">))</span> <span class="k">for</span> <span class="n">m</span> <span class="ow">in</span> <span class="n">matches</span><span class="p">)</span><span class="o">.</span><span class="n">types</span>
    <span class="n">left_fun</span> <span class="o">=</span> <span class="n">_KL_REGISTRY</span><span class="p">[</span><span class="n">left_p</span><span class="p">,</span> <span class="n">left_q</span><span class="p">]</span>
    <span class="n">right_fun</span> <span class="o">=</span> <span class="n">_KL_REGISTRY</span><span class="p">[</span><span class="n">right_p</span><span class="p">,</span> <span class="n">right_q</span><span class="p">]</span>
    <span class="k">if</span> <span class="n">left_fun</span> <span class="ow">is</span> <span class="ow">not</span> <span class="n">right_fun</span><span class="p">:</span>
        <span class="n">warnings</span><span class="o">.</span><span class="n">warn</span><span class="p">(</span><span class="s1">&#39;Ambiguous kl_divergence(</span><span class="si">{}</span><span class="s1">, </span><span class="si">{}</span><span class="s1">). Please register_kl(</span><span class="si">{}</span><span class="s1">, </span><span class="si">{}</span><span class="s1">)&#39;</span><span class="o">.</span><span class="n">format</span><span class="p">(</span>
            <span class="n">type_p</span><span class="o">.</span><span class="vm">__name__</span><span class="p">,</span> <span class="n">type_q</span><span class="o">.</span><span class="vm">__name__</span><span class="p">,</span> <span class="n">left_p</span><span class="o">.</span><span class="vm">__name__</span><span class="p">,</span> <span class="n">right_q</span><span class="o">.</span><span class="vm">__name__</span><span class="p">),</span>
            <span class="ne">RuntimeWarning</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">left_fun</span>


<span class="k">def</span> <span class="nf">_infinite_like</span><span class="p">(</span><span class="n">tensor</span><span class="p">):</span>
    <span class="sd">&quot;&quot;&quot;</span>
<span class="sd">    Helper function for obtaining infinite KL Divergence throughout</span>
<span class="sd">    &quot;&quot;&quot;</span>
    <span class="k">return</span> <span class="n">torch</span><span class="o">.</span><span class="n">full_like</span><span class="p">(</span><span class="n">tensor</span><span class="p">,</span> <span class="n">inf</span><span class="p">)</span>


<span class="k">def</span> <span class="nf">_x_log_x</span><span class="p">(</span><span class="n">tensor</span><span class="p">):</span>
    <span class="sd">&quot;&quot;&quot;</span>
<span class="sd">    Utility function for calculating x log x</span>
<span class="sd">    &quot;&quot;&quot;</span>
    <span class="k">return</span> <span class="n">tensor</span> <span class="o">*</span> <span class="n">tensor</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>


<span class="k">def</span> <span class="nf">_batch_trace_XXT</span><span class="p">(</span><span class="n">bmat</span><span class="p">):</span>
    <span class="sd">&quot;&quot;&quot;</span>
<span class="sd">    Utility function for calculating the trace of XX^{T} with X having arbitrary trailing batch dimensions</span>
<span class="sd">    &quot;&quot;&quot;</span>
    <span class="n">n</span> <span class="o">=</span> <span class="n">bmat</span><span class="o">.</span><span class="n">size</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
    <span class="n">m</span> <span class="o">=</span> <span class="n">bmat</span><span class="o">.</span><span class="n">size</span><span class="p">(</span><span class="o">-</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">flat_trace</span> <span class="o">=</span> <span class="n">bmat</span><span class="o">.</span><span class="n">reshape</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">,</span> <span class="n">m</span> <span class="o">*</span> <span class="n">n</span><span class="p">)</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">flat_trace</span><span class="o">.</span><span class="n">reshape</span><span class="p">(</span><span class="n">bmat</span><span class="o">.</span><span class="n">shape</span><span class="p">[:</span><span class="o">-</span><span class="mi">2</span><span class="p">])</span>


<div class="viewcode-block" id="kl_divergence"><a class="viewcode-back" href="../../../distributions.html#torch.distributions.kl.kl_divergence">[docs]</a><span class="k">def</span> <span class="nf">kl_divergence</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="sa">r</span><span class="sd">&quot;&quot;&quot;</span>
<span class="sd">    Compute Kullback-Leibler divergence :math:`KL(p \| q)` between two distributions.</span>

<span class="sd">    .. math::</span>

<span class="sd">        KL(p \| q) = \int p(x) \log\frac {p(x)} {q(x)} \,dx</span>

<span class="sd">    Args:</span>
<span class="sd">        p (Distribution): A :class:`~torch.distributions.Distribution` object.</span>
<span class="sd">        q (Distribution): A :class:`~torch.distributions.Distribution` object.</span>

<span class="sd">    Returns:</span>
<span class="sd">        Tensor: A batch of KL divergences of shape `batch_shape`.</span>

<span class="sd">    Raises:</span>
<span class="sd">        NotImplementedError: If the distribution types have not been registered via</span>
<span class="sd">            :meth:`register_kl`.</span>
<span class="sd">    &quot;&quot;&quot;</span>
    <span class="k">try</span><span class="p">:</span>
        <span class="n">fun</span> <span class="o">=</span> <span class="n">_KL_MEMOIZE</span><span class="p">[</span><span class="nb">type</span><span class="p">(</span><span class="n">p</span><span class="p">),</span> <span class="nb">type</span><span class="p">(</span><span class="n">q</span><span class="p">)]</span>
    <span class="k">except</span> <span class="ne">KeyError</span><span class="p">:</span>
        <span class="n">fun</span> <span class="o">=</span> <span class="n">_dispatch_kl</span><span class="p">(</span><span class="nb">type</span><span class="p">(</span><span class="n">p</span><span class="p">),</span> <span class="nb">type</span><span class="p">(</span><span class="n">q</span><span class="p">))</span>
        <span class="n">_KL_MEMOIZE</span><span class="p">[</span><span class="nb">type</span><span class="p">(</span><span class="n">p</span><span class="p">),</span> <span class="nb">type</span><span class="p">(</span><span class="n">q</span><span class="p">)]</span> <span class="o">=</span> <span class="n">fun</span>
    <span class="k">if</span> <span class="n">fun</span> <span class="ow">is</span> <span class="bp">NotImplemented</span><span class="p">:</span>
        <span class="k">raise</span> <span class="ne">NotImplementedError</span>
    <span class="k">return</span> <span class="n">fun</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">)</span></div>


<span class="c1">################################################################################</span>
<span class="c1"># KL Divergence Implementations</span>
<span class="c1">################################################################################</span>

<span class="n">_euler_gamma</span> <span class="o">=</span> <span class="mf">0.57721566490153286060</span>

<span class="c1"># Same distributions</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Bernoulli</span><span class="p">,</span> <span class="n">Bernoulli</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_bernoulli_bernoulli</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">probs</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">probs</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t1</span><span class="p">[</span><span class="n">q</span><span class="o">.</span><span class="n">probs</span> <span class="o">==</span> <span class="mi">0</span><span class="p">]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="n">t1</span><span class="p">[</span><span class="n">p</span><span class="o">.</span><span class="n">probs</span> <span class="o">==</span> <span class="mi">0</span><span class="p">]</span> <span class="o">=</span> <span class="mi">0</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="p">(</span><span class="mi">1</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span> <span class="o">*</span> <span class="p">((</span><span class="mi">1</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span> <span class="o">/</span> <span class="p">(</span><span class="mi">1</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">probs</span><span class="p">))</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span><span class="p">[</span><span class="n">q</span><span class="o">.</span><span class="n">probs</span> <span class="o">==</span> <span class="mi">1</span><span class="p">]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="n">t2</span><span class="p">[</span><span class="n">p</span><span class="o">.</span><span class="n">probs</span> <span class="o">==</span> <span class="mi">1</span><span class="p">]</span> <span class="o">=</span> <span class="mi">0</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Beta</span><span class="p">,</span> <span class="n">Beta</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_beta_beta</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">sum_params_p</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration0</span>
    <span class="n">sum_params_q</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">+</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration0</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration1</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">+</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration0</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">+</span> <span class="p">(</span><span class="n">sum_params_p</span><span class="p">)</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration1</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration0</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">+</span> <span class="p">(</span><span class="n">sum_params_q</span><span class="p">)</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration1</span><span class="p">)</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">digamma</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span><span class="p">)</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration0</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration0</span><span class="p">)</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">digamma</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration0</span><span class="p">)</span>
    <span class="n">t5</span> <span class="o">=</span> <span class="p">(</span><span class="n">sum_params_q</span> <span class="o">-</span> <span class="n">sum_params_p</span><span class="p">)</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">digamma</span><span class="p">(</span><span class="n">sum_params_p</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">-</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span> <span class="o">+</span> <span class="n">t5</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Binomial</span><span class="p">,</span> <span class="n">Binomial</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_binomial_binomial</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="c1"># from https://math.stackexchange.com/questions/2214993/</span>
    <span class="c1"># kullback-leibler-divergence-for-binomial-distributions-p-and-q</span>
    <span class="k">if</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">total_count</span> <span class="o">&lt;</span> <span class="n">q</span><span class="o">.</span><span class="n">total_count</span><span class="p">)</span><span class="o">.</span><span class="n">any</span><span class="p">():</span>
        <span class="k">raise</span> <span class="ne">NotImplementedError</span><span class="p">(</span><span class="s1">&#39;KL between Binomials where q.total_count &gt; p.total_count is not implemented&#39;</span><span class="p">)</span>
    <span class="n">kl</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">total_count</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">probs</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">logits</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">logits</span><span class="p">)</span> <span class="o">+</span> <span class="p">(</span><span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span><span class="o">.</span><span class="n">log1p</span><span class="p">()</span> <span class="o">-</span> <span class="p">(</span><span class="o">-</span><span class="n">q</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span><span class="o">.</span><span class="n">log1p</span><span class="p">())</span>
    <span class="n">inf_idxs</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">total_count</span> <span class="o">&gt;</span> <span class="n">q</span><span class="o">.</span><span class="n">total_count</span>
    <span class="n">kl</span><span class="p">[</span><span class="n">inf_idxs</span><span class="p">]</span> <span class="o">=</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">kl</span><span class="p">[</span><span class="n">inf_idxs</span><span class="p">])</span>
    <span class="k">return</span> <span class="n">kl</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Categorical</span><span class="p">,</span> <span class="n">Categorical</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_categorical_categorical</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">t</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">probs</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">logits</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">logits</span><span class="p">)</span>
    <span class="n">t</span><span class="p">[(</span><span class="n">q</span><span class="o">.</span><span class="n">probs</span> <span class="o">==</span> <span class="mi">0</span><span class="p">)</span><span class="o">.</span><span class="n">expand_as</span><span class="p">(</span><span class="n">t</span><span class="p">)]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="n">t</span><span class="p">[(</span><span class="n">p</span><span class="o">.</span><span class="n">probs</span> <span class="o">==</span> <span class="mi">0</span><span class="p">)</span><span class="o">.</span><span class="n">expand_as</span><span class="p">(</span><span class="n">t</span><span class="p">)]</span> <span class="o">=</span> <span class="mi">0</span>
    <span class="k">return</span> <span class="n">t</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">ContinuousBernoulli</span><span class="p">,</span> <span class="n">ContinuousBernoulli</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_continuous_bernoulli_continuous_bernoulli</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">mean</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">logits</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">logits</span><span class="p">)</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">_cont_bern_log_norm</span><span class="p">()</span> <span class="o">+</span> <span class="n">torch</span><span class="o">.</span><span class="n">log1p</span><span class="p">(</span><span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">_cont_bern_log_norm</span><span class="p">()</span> <span class="o">-</span> <span class="n">torch</span><span class="o">.</span><span class="n">log1p</span><span class="p">(</span><span class="o">-</span><span class="n">q</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Dirichlet</span><span class="p">,</span> <span class="n">Dirichlet</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_dirichlet_dirichlet</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="c1"># From http://bariskurt.com/kullback-leibler-divergence-between-two-dirichlet-and-beta-distributions/</span>
    <span class="n">sum_p_concentration</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
    <span class="n">sum_q_concentration</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">sum_p_concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">-</span> <span class="n">sum_q_concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">())</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">digamma</span><span class="p">()</span> <span class="o">-</span> <span class="n">sum_p_concentration</span><span class="o">.</span><span class="n">digamma</span><span class="p">()</span><span class="o">.</span><span class="n">unsqueeze</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">-</span> <span class="n">t2</span> <span class="o">+</span> <span class="p">(</span><span class="n">t3</span> <span class="o">*</span> <span class="n">t4</span><span class="p">)</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Exponential</span><span class="p">,</span> <span class="n">Exponential</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_exponential_exponential</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">rate_ratio</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="o">-</span><span class="n">rate_ratio</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">rate_ratio</span> <span class="o">-</span> <span class="mi">1</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">ExponentialFamily</span><span class="p">,</span> <span class="n">ExponentialFamily</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_expfamily_expfamily</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">if</span> <span class="ow">not</span> <span class="nb">type</span><span class="p">(</span><span class="n">p</span><span class="p">)</span> <span class="o">==</span> <span class="nb">type</span><span class="p">(</span><span class="n">q</span><span class="p">):</span>
        <span class="k">raise</span> <span class="ne">NotImplementedError</span><span class="p">(</span><span class="s2">&quot;The cross KL-divergence between different exponential families cannot </span><span class="se">\</span>
<span class="s2">                            be computed using Bregman divergences&quot;</span><span class="p">)</span>
    <span class="n">p_nparams</span> <span class="o">=</span> <span class="p">[</span><span class="n">np</span><span class="o">.</span><span class="n">detach</span><span class="p">()</span><span class="o">.</span><span class="n">requires_grad_</span><span class="p">()</span> <span class="k">for</span> <span class="n">np</span> <span class="ow">in</span> <span class="n">p</span><span class="o">.</span><span class="n">_natural_params</span><span class="p">]</span>
    <span class="n">q_nparams</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">_natural_params</span>
    <span class="n">lg_normal</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">_log_normalizer</span><span class="p">(</span><span class="o">*</span><span class="n">p_nparams</span><span class="p">)</span>
    <span class="n">gradients</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">autograd</span><span class="o">.</span><span class="n">grad</span><span class="p">(</span><span class="n">lg_normal</span><span class="o">.</span><span class="n">sum</span><span class="p">(),</span> <span class="n">p_nparams</span><span class="p">,</span> <span class="n">create_graph</span><span class="o">=</span><span class="kc">True</span><span class="p">)</span>
    <span class="n">result</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">_log_normalizer</span><span class="p">(</span><span class="o">*</span><span class="n">q_nparams</span><span class="p">)</span> <span class="o">-</span> <span class="n">lg_normal</span>
    <span class="k">for</span> <span class="n">pnp</span><span class="p">,</span> <span class="n">qnp</span><span class="p">,</span> <span class="n">g</span> <span class="ow">in</span> <span class="nb">zip</span><span class="p">(</span><span class="n">p_nparams</span><span class="p">,</span> <span class="n">q_nparams</span><span class="p">,</span> <span class="n">gradients</span><span class="p">):</span>
        <span class="n">term</span> <span class="o">=</span> <span class="p">(</span><span class="n">qnp</span> <span class="o">-</span> <span class="n">pnp</span><span class="p">)</span> <span class="o">*</span> <span class="n">g</span>
        <span class="n">result</span> <span class="o">-=</span> <span class="n">_sum_rightmost</span><span class="p">(</span><span class="n">term</span><span class="p">,</span> <span class="nb">len</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">event_shape</span><span class="p">))</span>
    <span class="k">return</span> <span class="n">result</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gamma</span><span class="p">,</span> <span class="n">Gamma</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_gamma_gamma</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">rate</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">lgamma</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="p">)</span> <span class="o">-</span> <span class="n">torch</span><span class="o">.</span><span class="n">lgamma</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="p">)</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="p">)</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">digamma</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="p">)</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gumbel</span><span class="p">,</span> <span class="n">Gumbel</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_gumbel_gumbel</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">ct1</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">ct2</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">ct3</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">loc</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="o">-</span><span class="n">ct1</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">-</span> <span class="n">ct2</span> <span class="o">+</span> <span class="n">ct3</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">ct1</span> <span class="o">*</span> <span class="n">_euler_gamma</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">exp</span><span class="p">(</span><span class="n">ct2</span> <span class="o">+</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="n">ct1</span><span class="p">)</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">-</span> <span class="n">ct3</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span> <span class="o">-</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="n">_euler_gamma</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Geometric</span><span class="p">,</span> <span class="n">Geometric</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_geometric_geometric</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span> <span class="o">-</span> <span class="n">torch</span><span class="o">.</span><span class="n">log1p</span><span class="p">(</span><span class="o">-</span><span class="n">q</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">probs</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">logits</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">HalfNormal</span><span class="p">,</span> <span class="n">HalfNormal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_halfnormal_halfnormal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_kl_normal_normal</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">base_dist</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">base_dist</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Laplace</span><span class="p">,</span> <span class="n">Laplace</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_laplace_laplace</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="c1"># From http://www.mast.queensu.ca/~communications/Papers/gil-msc11.pdf</span>
    <span class="n">scale_ratio</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">loc_abs_diff</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">loc</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span><span class="o">.</span><span class="n">abs</span><span class="p">()</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="o">-</span><span class="n">scale_ratio</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">loc_abs_diff</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">scale_ratio</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">exp</span><span class="p">(</span><span class="o">-</span><span class="n">loc_abs_diff</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span> <span class="o">-</span> <span class="mi">1</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">LowRankMultivariateNormal</span><span class="p">,</span> <span class="n">LowRankMultivariateNormal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_lowrankmultivariatenormal_lowrankmultivariatenormal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">if</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span> <span class="o">!=</span> <span class="n">q</span><span class="o">.</span><span class="n">event_shape</span><span class="p">:</span>
        <span class="k">raise</span> <span class="ne">ValueError</span><span class="p">(</span><span class="s2">&quot;KL-divergence between two Low Rank Multivariate Normals with</span><span class="se">\</span>
<span class="s2">                          different event shapes cannot be computed&quot;</span><span class="p">)</span>

    <span class="n">term1</span> <span class="o">=</span> <span class="p">(</span><span class="n">_batch_lowrank_logdet</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="p">,</span>
                                   <span class="n">q</span><span class="o">.</span><span class="n">_capacitance_tril</span><span class="p">)</span> <span class="o">-</span>
             <span class="n">_batch_lowrank_logdet</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="p">,</span> <span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="p">,</span>
                                   <span class="n">p</span><span class="o">.</span><span class="n">_capacitance_tril</span><span class="p">))</span>
    <span class="n">term3</span> <span class="o">=</span> <span class="n">_batch_lowrank_mahalanobis</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="p">,</span>
                                       <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">loc</span><span class="p">,</span>
                                       <span class="n">q</span><span class="o">.</span><span class="n">_capacitance_tril</span><span class="p">)</span>
    <span class="c1"># Expands term2 according to</span>
    <span class="c1"># inv(qcov) @ pcov = [inv(qD) - inv(qD) @ qW @ inv(qC) @ qW.T @ inv(qD)] @ (pW @ pW.T + pD)</span>
    <span class="c1">#                  = [inv(qD) - A.T @ A] @ (pD + pW @ pW.T)</span>
    <span class="n">qWt_qDinv</span> <span class="o">=</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="o">.</span><span class="n">transpose</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">,</span> <span class="o">-</span><span class="mi">2</span><span class="p">)</span> <span class="o">/</span>
                 <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="o">.</span><span class="n">unsqueeze</span><span class="p">(</span><span class="o">-</span><span class="mi">2</span><span class="p">))</span>
    <span class="n">A</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">triangular_solve</span><span class="p">(</span><span class="n">qWt_qDinv</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">_capacitance_tril</span><span class="p">,</span> <span class="n">upper</span><span class="o">=</span><span class="kc">False</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
    <span class="n">term21</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="p">)</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span>
    <span class="n">term22</span> <span class="o">=</span> <span class="n">_batch_trace_XXT</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span> <span class="o">*</span>
                              <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="o">.</span><span class="n">rsqrt</span><span class="p">()</span><span class="o">.</span><span class="n">unsqueeze</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">))</span>
    <span class="n">term23</span> <span class="o">=</span> <span class="n">_batch_trace_XXT</span><span class="p">(</span><span class="n">A</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="o">.</span><span class="n">sqrt</span><span class="p">()</span><span class="o">.</span><span class="n">unsqueeze</span><span class="p">(</span><span class="o">-</span><span class="mi">2</span><span class="p">))</span>
    <span class="n">term24</span> <span class="o">=</span> <span class="n">_batch_trace_XXT</span><span class="p">(</span><span class="n">A</span><span class="o">.</span><span class="n">matmul</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="p">))</span>
    <span class="n">term2</span> <span class="o">=</span> <span class="n">term21</span> <span class="o">+</span> <span class="n">term22</span> <span class="o">-</span> <span class="n">term23</span> <span class="o">-</span> <span class="n">term24</span>
    <span class="k">return</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">term1</span> <span class="o">+</span> <span class="n">term2</span> <span class="o">+</span> <span class="n">term3</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span><span class="p">[</span><span class="mi">0</span><span class="p">])</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">MultivariateNormal</span><span class="p">,</span> <span class="n">LowRankMultivariateNormal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_multivariatenormal_lowrankmultivariatenormal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">if</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span> <span class="o">!=</span> <span class="n">q</span><span class="o">.</span><span class="n">event_shape</span><span class="p">:</span>
        <span class="k">raise</span> <span class="ne">ValueError</span><span class="p">(</span><span class="s2">&quot;KL-divergence between two (Low Rank) Multivariate Normals with</span><span class="se">\</span>
<span class="s2">                          different event shapes cannot be computed&quot;</span><span class="p">)</span>

    <span class="n">term1</span> <span class="o">=</span> <span class="p">(</span><span class="n">_batch_lowrank_logdet</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="p">,</span>
                                   <span class="n">q</span><span class="o">.</span><span class="n">_capacitance_tril</span><span class="p">)</span> <span class="o">-</span>
             <span class="mi">2</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">diagonal</span><span class="p">(</span><span class="n">dim1</span><span class="o">=-</span><span class="mi">2</span><span class="p">,</span> <span class="n">dim2</span><span class="o">=-</span><span class="mi">1</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">))</span>
    <span class="n">term3</span> <span class="o">=</span> <span class="n">_batch_lowrank_mahalanobis</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="p">,</span>
                                       <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">loc</span><span class="p">,</span>
                                       <span class="n">q</span><span class="o">.</span><span class="n">_capacitance_tril</span><span class="p">)</span>
    <span class="c1"># Expands term2 according to</span>
    <span class="c1"># inv(qcov) @ pcov = [inv(qD) - inv(qD) @ qW @ inv(qC) @ qW.T @ inv(qD)] @ p_tril @ p_tril.T</span>
    <span class="c1">#                  = [inv(qD) - A.T @ A] @ p_tril @ p_tril.T</span>
    <span class="n">qWt_qDinv</span> <span class="o">=</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="o">.</span><span class="n">transpose</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">,</span> <span class="o">-</span><span class="mi">2</span><span class="p">)</span> <span class="o">/</span>
                 <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="o">.</span><span class="n">unsqueeze</span><span class="p">(</span><span class="o">-</span><span class="mi">2</span><span class="p">))</span>
    <span class="n">A</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">triangular_solve</span><span class="p">(</span><span class="n">qWt_qDinv</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">_capacitance_tril</span><span class="p">,</span> <span class="n">upper</span><span class="o">=</span><span class="kc">False</span><span class="p">)[</span><span class="mi">0</span><span class="p">]</span>
    <span class="n">term21</span> <span class="o">=</span> <span class="n">_batch_trace_XXT</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span> <span class="o">*</span>
                              <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="o">.</span><span class="n">rsqrt</span><span class="p">()</span><span class="o">.</span><span class="n">unsqueeze</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">))</span>
    <span class="n">term22</span> <span class="o">=</span> <span class="n">_batch_trace_XXT</span><span class="p">(</span><span class="n">A</span><span class="o">.</span><span class="n">matmul</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="p">))</span>
    <span class="n">term2</span> <span class="o">=</span> <span class="n">term21</span> <span class="o">-</span> <span class="n">term22</span>
    <span class="k">return</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">term1</span> <span class="o">+</span> <span class="n">term2</span> <span class="o">+</span> <span class="n">term3</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span><span class="p">[</span><span class="mi">0</span><span class="p">])</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">LowRankMultivariateNormal</span><span class="p">,</span> <span class="n">MultivariateNormal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_lowrankmultivariatenormal_multivariatenormal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">if</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span> <span class="o">!=</span> <span class="n">q</span><span class="o">.</span><span class="n">event_shape</span><span class="p">:</span>
        <span class="k">raise</span> <span class="ne">ValueError</span><span class="p">(</span><span class="s2">&quot;KL-divergence between two (Low Rank) Multivariate Normals with</span><span class="se">\</span>
<span class="s2">                          different event shapes cannot be computed&quot;</span><span class="p">)</span>

    <span class="n">term1</span> <span class="o">=</span> <span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">diagonal</span><span class="p">(</span><span class="n">dim1</span><span class="o">=-</span><span class="mi">2</span><span class="p">,</span> <span class="n">dim2</span><span class="o">=-</span><span class="mi">1</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span> <span class="o">-</span>
             <span class="n">_batch_lowrank_logdet</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="p">,</span> <span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="p">,</span>
                                   <span class="n">p</span><span class="o">.</span><span class="n">_capacitance_tril</span><span class="p">))</span>
    <span class="n">term3</span> <span class="o">=</span> <span class="n">_batch_mahalanobis</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="p">,</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">loc</span><span class="p">))</span>
    <span class="c1"># Expands term2 according to</span>
    <span class="c1"># inv(qcov) @ pcov = inv(q_tril @ q_tril.T) @ (pW @ pW.T + pD)</span>
    <span class="n">combined_batch_shape</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">_C</span><span class="o">.</span><span class="n">_infer_size</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">shape</span><span class="p">[:</span><span class="o">-</span><span class="mi">2</span><span class="p">],</span>
                                                <span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="o">.</span><span class="n">shape</span><span class="p">[:</span><span class="o">-</span><span class="mi">2</span><span class="p">])</span>
    <span class="n">n</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
    <span class="n">q_scale_tril</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">expand</span><span class="p">(</span><span class="n">combined_batch_shape</span> <span class="o">+</span> <span class="p">(</span><span class="n">n</span><span class="p">,</span> <span class="n">n</span><span class="p">))</span>
    <span class="n">p_cov_factor</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_factor</span><span class="o">.</span><span class="n">expand</span><span class="p">(</span><span class="n">combined_batch_shape</span> <span class="o">+</span>
                                                      <span class="p">(</span><span class="n">n</span><span class="p">,</span> <span class="n">p</span><span class="o">.</span><span class="n">cov_factor</span><span class="o">.</span><span class="n">size</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)))</span>
    <span class="n">p_cov_diag</span> <span class="o">=</span> <span class="p">(</span><span class="n">torch</span><span class="o">.</span><span class="n">diag_embed</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_cov_diag</span><span class="o">.</span><span class="n">sqrt</span><span class="p">())</span>
                  <span class="o">.</span><span class="n">expand</span><span class="p">(</span><span class="n">combined_batch_shape</span> <span class="o">+</span> <span class="p">(</span><span class="n">n</span><span class="p">,</span> <span class="n">n</span><span class="p">)))</span>
    <span class="n">term21</span> <span class="o">=</span> <span class="n">_batch_trace_XXT</span><span class="p">(</span><span class="n">torch</span><span class="o">.</span><span class="n">triangular_solve</span><span class="p">(</span><span class="n">p_cov_factor</span><span class="p">,</span> <span class="n">q_scale_tril</span><span class="p">,</span> <span class="n">upper</span><span class="o">=</span><span class="kc">False</span><span class="p">)[</span><span class="mi">0</span><span class="p">])</span>
    <span class="n">term22</span> <span class="o">=</span> <span class="n">_batch_trace_XXT</span><span class="p">(</span><span class="n">torch</span><span class="o">.</span><span class="n">triangular_solve</span><span class="p">(</span><span class="n">p_cov_diag</span><span class="p">,</span> <span class="n">q_scale_tril</span><span class="p">,</span> <span class="n">upper</span><span class="o">=</span><span class="kc">False</span><span class="p">)[</span><span class="mi">0</span><span class="p">])</span>
    <span class="n">term2</span> <span class="o">=</span> <span class="n">term21</span> <span class="o">+</span> <span class="n">term22</span>
    <span class="k">return</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">term1</span> <span class="o">+</span> <span class="n">term2</span> <span class="o">+</span> <span class="n">term3</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span><span class="p">[</span><span class="mi">0</span><span class="p">])</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">MultivariateNormal</span><span class="p">,</span> <span class="n">MultivariateNormal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_multivariatenormal_multivariatenormal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="c1"># From https://en.wikipedia.org/wiki/Multivariate_normal_distribution#Kullback%E2%80%93Leibler_divergence</span>
    <span class="k">if</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span> <span class="o">!=</span> <span class="n">q</span><span class="o">.</span><span class="n">event_shape</span><span class="p">:</span>
        <span class="k">raise</span> <span class="ne">ValueError</span><span class="p">(</span><span class="s2">&quot;KL-divergence between two Multivariate Normals with</span><span class="se">\</span>
<span class="s2">                          different event shapes cannot be computed&quot;</span><span class="p">)</span>

    <span class="n">half_term1</span> <span class="o">=</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">diagonal</span><span class="p">(</span><span class="n">dim1</span><span class="o">=-</span><span class="mi">2</span><span class="p">,</span> <span class="n">dim2</span><span class="o">=-</span><span class="mi">1</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">)</span> <span class="o">-</span>
                  <span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">diagonal</span><span class="p">(</span><span class="n">dim1</span><span class="o">=-</span><span class="mi">2</span><span class="p">,</span> <span class="n">dim2</span><span class="o">=-</span><span class="mi">1</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span><span class="o">.</span><span class="n">sum</span><span class="p">(</span><span class="o">-</span><span class="mi">1</span><span class="p">))</span>
    <span class="n">combined_batch_shape</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">_C</span><span class="o">.</span><span class="n">_infer_size</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">shape</span><span class="p">[:</span><span class="o">-</span><span class="mi">2</span><span class="p">],</span>
                                                <span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">shape</span><span class="p">[:</span><span class="o">-</span><span class="mi">2</span><span class="p">])</span>
    <span class="n">n</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span><span class="p">[</span><span class="mi">0</span><span class="p">]</span>
    <span class="n">q_scale_tril</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">expand</span><span class="p">(</span><span class="n">combined_batch_shape</span> <span class="o">+</span> <span class="p">(</span><span class="n">n</span><span class="p">,</span> <span class="n">n</span><span class="p">))</span>
    <span class="n">p_scale_tril</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="o">.</span><span class="n">expand</span><span class="p">(</span><span class="n">combined_batch_shape</span> <span class="o">+</span> <span class="p">(</span><span class="n">n</span><span class="p">,</span> <span class="n">n</span><span class="p">))</span>
    <span class="n">term2</span> <span class="o">=</span> <span class="n">_batch_trace_XXT</span><span class="p">(</span><span class="n">torch</span><span class="o">.</span><span class="n">triangular_solve</span><span class="p">(</span><span class="n">p_scale_tril</span><span class="p">,</span> <span class="n">q_scale_tril</span><span class="p">,</span> <span class="n">upper</span><span class="o">=</span><span class="kc">False</span><span class="p">)[</span><span class="mi">0</span><span class="p">])</span>
    <span class="n">term3</span> <span class="o">=</span> <span class="n">_batch_mahalanobis</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">_unbroadcasted_scale_tril</span><span class="p">,</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">loc</span><span class="p">))</span>
    <span class="k">return</span> <span class="n">half_term1</span> <span class="o">+</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">term2</span> <span class="o">+</span> <span class="n">term3</span> <span class="o">-</span> <span class="n">n</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Normal</span><span class="p">,</span> <span class="n">Normal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_normal_normal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">var_ratio</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="p">)</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="p">((</span><span class="n">p</span><span class="o">.</span><span class="n">loc</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="p">)</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="k">return</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">var_ratio</span> <span class="o">+</span> <span class="n">t1</span> <span class="o">-</span> <span class="mi">1</span> <span class="o">-</span> <span class="n">var_ratio</span><span class="o">.</span><span class="n">log</span><span class="p">())</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">OneHotCategorical</span><span class="p">,</span> <span class="n">OneHotCategorical</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_onehotcategorical_onehotcategorical</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_kl_categorical_categorical</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">_categorical</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">_categorical</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Pareto</span><span class="p">,</span> <span class="n">Pareto</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_pareto_pareto</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="c1"># From http://www.mast.queensu.ca/~communications/Papers/gil-msc11.pdf</span>
    <span class="n">scale_ratio</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">alpha_ratio</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">alpha</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">alpha</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">alpha</span> <span class="o">*</span> <span class="n">scale_ratio</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="o">-</span><span class="n">alpha_ratio</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">result</span> <span class="o">=</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">alpha_ratio</span> <span class="o">-</span> <span class="mi">1</span>
    <span class="n">result</span><span class="p">[</span><span class="n">p</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">lower_bound</span> <span class="o">&lt;</span> <span class="n">q</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">lower_bound</span><span class="p">]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Poisson</span><span class="p">,</span> <span class="n">Poisson</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_poisson_poisson</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">log</span><span class="p">())</span> <span class="o">-</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">rate</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">TransformedDistribution</span><span class="p">,</span> <span class="n">TransformedDistribution</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_transformed_transformed</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">if</span> <span class="n">p</span><span class="o">.</span><span class="n">transforms</span> <span class="o">!=</span> <span class="n">q</span><span class="o">.</span><span class="n">transforms</span><span class="p">:</span>
        <span class="k">raise</span> <span class="ne">NotImplementedError</span>
    <span class="k">if</span> <span class="n">p</span><span class="o">.</span><span class="n">event_shape</span> <span class="o">!=</span> <span class="n">q</span><span class="o">.</span><span class="n">event_shape</span><span class="p">:</span>
        <span class="k">raise</span> <span class="ne">NotImplementedError</span>
    <span class="c1"># extra_event_dim = len(p.event_shape) - len(p.base_dist.event_shape)</span>
    <span class="n">extra_event_dim</span> <span class="o">=</span> <span class="nb">len</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">event_shape</span><span class="p">)</span>
    <span class="n">base_kl_divergence</span> <span class="o">=</span> <span class="n">kl_divergence</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">base_dist</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">base_dist</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">_sum_rightmost</span><span class="p">(</span><span class="n">base_kl_divergence</span><span class="p">,</span> <span class="n">extra_event_dim</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Uniform</span><span class="p">,</span> <span class="n">Uniform</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_uniform_uniform</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">result</span> <span class="o">=</span> <span class="p">((</span><span class="n">q</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">low</span><span class="p">)</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">))</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">result</span><span class="p">[(</span><span class="n">q</span><span class="o">.</span><span class="n">low</span> <span class="o">&gt;</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">)</span> <span class="o">|</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">high</span> <span class="o">&lt;</span> <span class="n">p</span><span class="o">.</span><span class="n">high</span><span class="p">)]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>


<span class="c1"># Different distributions</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Bernoulli</span><span class="p">,</span> <span class="n">Poisson</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_bernoulli_poisson</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span> <span class="o">-</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">probs</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Beta</span><span class="p">,</span> <span class="n">ContinuousBernoulli</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_beta_continuous_bernoulli</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">mean</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">logits</span> <span class="o">-</span> <span class="n">torch</span><span class="o">.</span><span class="n">log1p</span><span class="p">(</span><span class="o">-</span><span class="n">q</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">_cont_bern_log_norm</span><span class="p">()</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Beta</span><span class="p">,</span> <span class="n">Pareto</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_beta_infinity</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Beta</span><span class="p">,</span> <span class="n">Exponential</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_beta_exponential</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">+</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration0</span><span class="p">))</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Beta</span><span class="p">,</span> <span class="n">Gamma</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_beta_gamma</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">concentration</span> <span class="o">-</span> <span class="mi">1</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span><span class="o">.</span><span class="n">digamma</span><span class="p">()</span> <span class="o">-</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration0</span><span class="p">)</span><span class="o">.</span><span class="n">digamma</span><span class="p">())</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration0</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">-</span> <span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span>

<span class="c1"># TODO: Add Beta-Laplace KL Divergence</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Beta</span><span class="p">,</span> <span class="n">Normal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_beta_normal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">E_beta</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration0</span><span class="p">)</span>
    <span class="n">var_normal</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">var_normal</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">(</span><span class="n">E_beta</span> <span class="o">*</span> <span class="p">(</span><span class="mi">1</span> <span class="o">-</span> <span class="n">E_beta</span><span class="p">)</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration0</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span> <span class="o">+</span> <span class="n">E_beta</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">))</span> <span class="o">*</span> <span class="mf">0.5</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">*</span> <span class="n">E_beta</span>
    <span class="n">t5</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span> <span class="o">*</span> <span class="mf">0.5</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="p">(</span><span class="n">t3</span> <span class="o">-</span> <span class="n">t4</span> <span class="o">+</span> <span class="n">t5</span><span class="p">)</span> <span class="o">/</span> <span class="n">var_normal</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Beta</span><span class="p">,</span> <span class="n">Uniform</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_beta_uniform</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">result</span> <span class="o">=</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span> <span class="o">+</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">low</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">result</span><span class="p">[(</span><span class="n">q</span><span class="o">.</span><span class="n">low</span> <span class="o">&gt;</span> <span class="n">p</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">lower_bound</span><span class="p">)</span> <span class="o">|</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">high</span> <span class="o">&lt;</span> <span class="n">p</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">upper_bound</span><span class="p">)]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>

<span class="c1"># Note that the KL between a ContinuousBernoulli and Beta has no closed form</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">ContinuousBernoulli</span><span class="p">,</span> <span class="n">Pareto</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_continuous_bernoulli_infinity</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">ContinuousBernoulli</span><span class="p">,</span> <span class="n">Exponential</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_continuous_bernoulli_exponential</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span> <span class="o">-</span> <span class="n">torch</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="p">)</span> <span class="o">+</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">mean</span>

<span class="c1"># Note that the KL between a ContinuousBernoulli and Gamma has no closed form</span>
<span class="c1"># TODO: Add ContinuousBernoulli-Laplace KL Divergence</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">ContinuousBernoulli</span><span class="p">,</span> <span class="n">Normal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_continuous_bernoulli_normal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">math</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="mf">2.</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span><span class="p">)</span> <span class="o">+</span> <span class="n">torch</span><span class="o">.</span><span class="n">square</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="p">))</span> <span class="o">+</span> <span class="n">torch</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="p">)</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">variance</span> <span class="o">+</span> <span class="n">torch</span><span class="o">.</span><span class="n">square</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">mean</span><span class="p">)</span> <span class="o">-</span> <span class="mf">2.</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">mean</span><span class="p">)</span> <span class="o">/</span> <span class="p">(</span><span class="mf">2.0</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">square</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="p">))</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">ContinuousBernoulli</span><span class="p">,</span> <span class="n">Uniform</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_continuous_bernoulli_uniform</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">result</span> <span class="o">=</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span> <span class="o">+</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">low</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="k">return</span> <span class="n">torch</span><span class="o">.</span><span class="n">where</span><span class="p">(</span><span class="n">torch</span><span class="o">.</span><span class="n">max</span><span class="p">(</span><span class="n">torch</span><span class="o">.</span><span class="n">ge</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">low</span><span class="p">,</span> <span class="n">p</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">lower_bound</span><span class="p">),</span>
                                 <span class="n">torch</span><span class="o">.</span><span class="n">le</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">high</span><span class="p">,</span> <span class="n">p</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">upper_bound</span><span class="p">)),</span>
                       <span class="n">torch</span><span class="o">.</span><span class="n">ones_like</span><span class="p">(</span><span class="n">result</span><span class="p">)</span> <span class="o">*</span> <span class="n">inf</span><span class="p">,</span> <span class="n">result</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Exponential</span><span class="p">,</span> <span class="n">Beta</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Exponential</span><span class="p">,</span> <span class="n">ContinuousBernoulli</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Exponential</span><span class="p">,</span> <span class="n">Pareto</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Exponential</span><span class="p">,</span> <span class="n">Uniform</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_exponential_infinity</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">rate</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Exponential</span><span class="p">,</span> <span class="n">Gamma</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_exponential_gamma</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">ratio</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="o">-</span><span class="n">q</span><span class="o">.</span><span class="n">concentration</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="n">ratio</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">ratio</span> <span class="o">+</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">+</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span> <span class="o">*</span> <span class="n">_euler_gamma</span> <span class="o">-</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="n">_euler_gamma</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Exponential</span><span class="p">,</span> <span class="n">Gumbel</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_exponential_gumbel</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">scale_rate_prod</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">loc_scale_ratio</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">scale_rate_prod</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">-</span> <span class="mi">1</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">exp</span><span class="p">(</span><span class="n">loc_scale_ratio</span><span class="p">)</span> <span class="o">*</span> <span class="n">scale_rate_prod</span> <span class="o">/</span> <span class="p">(</span><span class="n">scale_rate_prod</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">scale_rate_prod</span><span class="o">.</span><span class="n">reciprocal</span><span class="p">()</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">-</span> <span class="n">loc_scale_ratio</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span>

<span class="c1"># TODO: Add Exponential-Laplace KL Divergence</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Exponential</span><span class="p">,</span> <span class="n">Normal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_exponential_normal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">var_normal</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">rate_sqr</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="n">rate_sqr</span> <span class="o">*</span> <span class="n">var_normal</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span><span class="p">)</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">rate_sqr</span><span class="o">.</span><span class="n">reciprocal</span><span class="p">()</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span> <span class="o">*</span> <span class="mf">0.5</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">-</span> <span class="mi">1</span> <span class="o">+</span> <span class="p">(</span><span class="n">t2</span> <span class="o">-</span> <span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span><span class="p">)</span> <span class="o">/</span> <span class="n">var_normal</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gamma</span><span class="p">,</span> <span class="n">Beta</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gamma</span><span class="p">,</span> <span class="n">ContinuousBernoulli</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gamma</span><span class="p">,</span> <span class="n">Pareto</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gamma</span><span class="p">,</span> <span class="n">Uniform</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_gamma_infinity</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gamma</span><span class="p">,</span> <span class="n">Exponential</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_gamma_exponential</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">+</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gamma</span><span class="p">,</span> <span class="n">Gumbel</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_gamma_gumbel</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">beta_scale_prod</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">loc_scale_ratio</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span> <span class="o">-</span> <span class="mi">1</span><span class="p">)</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">digamma</span><span class="p">()</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">beta_scale_prod</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span> <span class="o">/</span> <span class="n">beta_scale_prod</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">exp</span><span class="p">(</span><span class="n">loc_scale_ratio</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="n">beta_scale_prod</span><span class="o">.</span><span class="n">reciprocal</span><span class="p">())</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="p">)</span> <span class="o">-</span> <span class="n">loc_scale_ratio</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span>

<span class="c1"># TODO: Add Gamma-Laplace KL Divergence</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gamma</span><span class="p">,</span> <span class="n">Normal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_gamma_normal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">var_normal</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">beta_sqr</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="n">beta_sqr</span> <span class="o">*</span> <span class="n">var_normal</span> <span class="o">*</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span><span class="p">)</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="p">)</span> <span class="o">/</span> <span class="n">beta_sqr</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">rate</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">concentration</span> <span class="o">-</span> <span class="mi">1</span><span class="p">)</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">digamma</span><span class="p">()</span> <span class="o">+</span> <span class="p">(</span><span class="n">t2</span> <span class="o">-</span> <span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span><span class="p">)</span> <span class="o">/</span> <span class="n">var_normal</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gumbel</span><span class="p">,</span> <span class="n">Beta</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gumbel</span><span class="p">,</span> <span class="n">ContinuousBernoulli</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gumbel</span><span class="p">,</span> <span class="n">Exponential</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gumbel</span><span class="p">,</span> <span class="n">Gamma</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gumbel</span><span class="p">,</span> <span class="n">Pareto</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gumbel</span><span class="p">,</span> <span class="n">Uniform</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_gumbel_infinity</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span>

<span class="c1"># TODO: Add Gumbel-Laplace KL Divergence</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Gumbel</span><span class="p">,</span> <span class="n">Normal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_gumbel_normal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">param_ratio</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="p">(</span><span class="n">param_ratio</span> <span class="o">/</span> <span class="n">math</span><span class="o">.</span><span class="n">sqrt</span><span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span><span class="p">))</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="p">(</span><span class="n">math</span><span class="o">.</span><span class="n">pi</span> <span class="o">*</span> <span class="n">param_ratio</span> <span class="o">*</span> <span class="mf">0.5</span><span class="p">)</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span> <span class="o">/</span> <span class="mi">3</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">((</span><span class="n">p</span><span class="o">.</span><span class="n">loc</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">*</span> <span class="n">_euler_gamma</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="p">)</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span> <span class="o">*</span> <span class="mf">0.5</span>
    <span class="k">return</span> <span class="o">-</span><span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span> <span class="o">-</span> <span class="p">(</span><span class="n">_euler_gamma</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Laplace</span><span class="p">,</span> <span class="n">Beta</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Laplace</span><span class="p">,</span> <span class="n">ContinuousBernoulli</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Laplace</span><span class="p">,</span> <span class="n">Exponential</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Laplace</span><span class="p">,</span> <span class="n">Gamma</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Laplace</span><span class="p">,</span> <span class="n">Pareto</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Laplace</span><span class="p">,</span> <span class="n">Uniform</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_laplace_infinity</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Laplace</span><span class="p">,</span> <span class="n">Normal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_laplace_normal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">var_normal</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">scale_sqr_var_ratio</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span> <span class="o">/</span> <span class="n">var_normal</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="n">torch</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="n">scale_sqr_var_ratio</span> <span class="o">/</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span><span class="p">)</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">loc</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">loc</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="k">return</span> <span class="o">-</span><span class="n">t1</span> <span class="o">+</span> <span class="n">scale_sqr_var_ratio</span> <span class="o">+</span> <span class="p">(</span><span class="n">t2</span> <span class="o">-</span> <span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span><span class="p">)</span> <span class="o">/</span> <span class="n">var_normal</span> <span class="o">-</span> <span class="mi">1</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Normal</span><span class="p">,</span> <span class="n">Beta</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Normal</span><span class="p">,</span> <span class="n">ContinuousBernoulli</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Normal</span><span class="p">,</span> <span class="n">Exponential</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Normal</span><span class="p">,</span> <span class="n">Gamma</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Normal</span><span class="p">,</span> <span class="n">Pareto</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Normal</span><span class="p">,</span> <span class="n">Uniform</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_normal_infinity</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Normal</span><span class="p">,</span> <span class="n">Gumbel</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_normal_gumbel</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">mean_scale_ratio</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">loc</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">var_scale_sqr_ratio</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="p">)</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">loc_scale_ratio</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">var_scale_sqr_ratio</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">*</span> <span class="mf">0.5</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">mean_scale_ratio</span> <span class="o">-</span> <span class="n">loc_scale_ratio</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">exp</span><span class="p">(</span><span class="o">-</span><span class="n">mean_scale_ratio</span> <span class="o">+</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="n">var_scale_sqr_ratio</span> <span class="o">+</span> <span class="n">loc_scale_ratio</span><span class="p">)</span>
    <span class="k">return</span> <span class="o">-</span><span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span> <span class="o">-</span> <span class="p">(</span><span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="mi">1</span> <span class="o">+</span> <span class="n">math</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span><span class="p">)))</span>

<span class="c1"># TODO: Add Normal-Laplace KL Divergence</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Pareto</span><span class="p">,</span> <span class="n">Beta</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Pareto</span><span class="p">,</span> <span class="n">ContinuousBernoulli</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Pareto</span><span class="p">,</span> <span class="n">Uniform</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_pareto_infinity</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">scale</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Pareto</span><span class="p">,</span> <span class="n">Exponential</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_pareto_exponential</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">scale_rate_prod</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">/</span> <span class="n">scale_rate_prod</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">alpha</span><span class="o">.</span><span class="n">reciprocal</span><span class="p">()</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">*</span> <span class="n">scale_rate_prod</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">-</span> <span class="mi">1</span><span class="p">)</span>
    <span class="n">result</span> <span class="o">=</span> <span class="n">t1</span> <span class="o">-</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span> <span class="o">-</span> <span class="mi">1</span>
    <span class="n">result</span><span class="p">[</span><span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">&lt;=</span> <span class="mi">1</span><span class="p">]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Pareto</span><span class="p">,</span> <span class="n">Gamma</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_pareto_gamma</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">common_term</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">alpha</span><span class="o">.</span><span class="n">reciprocal</span><span class="p">()</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">alpha</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">-</span> <span class="n">common_term</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">(</span><span class="mi">1</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="p">)</span> <span class="o">*</span> <span class="n">common_term</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">-</span> <span class="mi">1</span><span class="p">)</span>
    <span class="n">result</span> <span class="o">=</span> <span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span> <span class="o">-</span> <span class="mi">1</span>
    <span class="n">result</span><span class="p">[</span><span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">&lt;=</span> <span class="mi">1</span><span class="p">]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>

<span class="c1"># TODO: Add Pareto-Laplace KL Divergence</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Pareto</span><span class="p">,</span> <span class="n">Normal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_pareto_normal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">var_normal</span> <span class="o">=</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">common_term</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">-</span> <span class="mi">1</span><span class="p">)</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="p">(</span><span class="n">math</span><span class="o">.</span><span class="n">sqrt</span><span class="p">(</span><span class="mi">2</span> <span class="o">*</span> <span class="n">math</span><span class="o">.</span><span class="n">pi</span><span class="p">)</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span> <span class="o">*</span> <span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">/</span> <span class="n">p</span><span class="o">.</span><span class="n">scale</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">alpha</span><span class="o">.</span><span class="n">reciprocal</span><span class="p">()</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">*</span> <span class="n">common_term</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">-</span> <span class="mi">2</span><span class="p">)</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">*</span> <span class="n">common_term</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="n">result</span> <span class="o">=</span> <span class="n">t1</span> <span class="o">-</span> <span class="n">t2</span> <span class="o">+</span> <span class="p">(</span><span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span><span class="p">)</span> <span class="o">/</span> <span class="n">var_normal</span> <span class="o">-</span> <span class="mi">1</span>
    <span class="n">result</span><span class="p">[</span><span class="n">p</span><span class="o">.</span><span class="n">alpha</span> <span class="o">&lt;=</span> <span class="mi">2</span><span class="p">]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Poisson</span><span class="p">,</span> <span class="n">Bernoulli</span><span class="p">)</span>
<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Poisson</span><span class="p">,</span> <span class="n">Binomial</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_poisson_infinity</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">return</span> <span class="n">_infinite_like</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">rate</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Uniform</span><span class="p">,</span> <span class="n">Beta</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_uniform_beta</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">common_term</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">torch</span><span class="o">.</span><span class="n">log</span><span class="p">(</span><span class="n">common_term</span><span class="p">)</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">-</span> <span class="mi">1</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="n">_x_log_x</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span><span class="p">)</span> <span class="o">-</span> <span class="n">_x_log_x</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">)</span> <span class="o">-</span> <span class="n">common_term</span><span class="p">)</span> <span class="o">/</span> <span class="n">common_term</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">concentration0</span> <span class="o">-</span> <span class="mi">1</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="n">_x_log_x</span><span class="p">((</span><span class="mi">1</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">high</span><span class="p">))</span> <span class="o">-</span> <span class="n">_x_log_x</span><span class="p">((</span><span class="mi">1</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">))</span> <span class="o">+</span> <span class="n">common_term</span><span class="p">)</span> <span class="o">/</span> <span class="n">common_term</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration1</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">+</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration0</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">-</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">concentration1</span> <span class="o">+</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration0</span><span class="p">)</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span>
    <span class="n">result</span> <span class="o">=</span> <span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span> <span class="o">-</span> <span class="n">t1</span> <span class="o">-</span> <span class="n">t2</span>
    <span class="n">result</span><span class="p">[(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">&gt;</span> <span class="n">q</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">upper_bound</span><span class="p">)</span> <span class="o">|</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">low</span> <span class="o">&lt;</span> <span class="n">q</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">lower_bound</span><span class="p">)]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Uniform</span><span class="p">,</span> <span class="n">ContinuousBernoulli</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_uniform_continuous_bernoulli</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">result</span> <span class="o">=</span> <span class="o">-</span><span class="n">p</span><span class="o">.</span><span class="n">entropy</span><span class="p">()</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">mean</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">logits</span> <span class="o">-</span> <span class="n">torch</span><span class="o">.</span><span class="n">log1p</span><span class="p">(</span><span class="o">-</span><span class="n">q</span><span class="o">.</span><span class="n">probs</span><span class="p">)</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">_cont_bern_log_norm</span><span class="p">()</span>
    <span class="k">return</span> <span class="n">torch</span><span class="o">.</span><span class="n">where</span><span class="p">(</span><span class="n">torch</span><span class="o">.</span><span class="n">max</span><span class="p">(</span><span class="n">torch</span><span class="o">.</span><span class="n">ge</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">upper_bound</span><span class="p">),</span>
                                 <span class="n">torch</span><span class="o">.</span><span class="n">le</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">lower_bound</span><span class="p">)),</span>
                       <span class="n">torch</span><span class="o">.</span><span class="n">ones_like</span><span class="p">(</span><span class="n">result</span><span class="p">)</span> <span class="o">*</span> <span class="n">inf</span><span class="p">,</span> <span class="n">result</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Uniform</span><span class="p">,</span> <span class="n">Exponential</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_uniform_exponetial</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">result</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span> <span class="o">-</span> <span class="p">((</span><span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">)</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">result</span><span class="p">[</span><span class="n">p</span><span class="o">.</span><span class="n">low</span> <span class="o">&lt;</span> <span class="n">q</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">lower_bound</span><span class="p">]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Uniform</span><span class="p">,</span> <span class="n">Gamma</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_uniform_gamma</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">common_term</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">common_term</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="o">.</span><span class="n">lgamma</span><span class="p">()</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">(</span><span class="mi">1</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">concentration</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="n">_x_log_x</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span><span class="p">)</span> <span class="o">-</span> <span class="n">_x_log_x</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">)</span> <span class="o">-</span> <span class="n">common_term</span><span class="p">)</span> <span class="o">/</span> <span class="n">common_term</span>
    <span class="n">t4</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">rate</span> <span class="o">*</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span>
    <span class="n">result</span> <span class="o">=</span> <span class="o">-</span><span class="n">t1</span> <span class="o">+</span> <span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span> <span class="o">+</span> <span class="n">t4</span>
    <span class="n">result</span><span class="p">[</span><span class="n">p</span><span class="o">.</span><span class="n">low</span> <span class="o">&lt;</span> <span class="n">q</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">lower_bound</span><span class="p">]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Uniform</span><span class="p">,</span> <span class="n">Gumbel</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_uniform_gumbel</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">common_term</span> <span class="o">=</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">)</span>
    <span class="n">high_loc_diff</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">low_loc_diff</span> <span class="o">=</span> <span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">low</span> <span class="o">-</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="n">common_term</span><span class="o">.</span><span class="n">log</span><span class="p">()</span> <span class="o">+</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">high_loc_diff</span> <span class="o">+</span> <span class="n">low_loc_diff</span><span class="p">)</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="n">common_term</span> <span class="o">*</span> <span class="p">(</span><span class="n">torch</span><span class="o">.</span><span class="n">exp</span><span class="p">(</span><span class="o">-</span><span class="n">high_loc_diff</span><span class="p">)</span> <span class="o">-</span> <span class="n">torch</span><span class="o">.</span><span class="n">exp</span><span class="p">(</span><span class="o">-</span><span class="n">low_loc_diff</span><span class="p">))</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">-</span> <span class="n">t2</span>

<span class="c1"># TODO: Uniform-Laplace KL Divergence</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Uniform</span><span class="p">,</span> <span class="n">Normal</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_uniform_normal</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">common_term</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="p">(</span><span class="n">math</span><span class="o">.</span><span class="n">sqrt</span><span class="p">(</span><span class="n">math</span><span class="o">.</span><span class="n">pi</span> <span class="o">*</span> <span class="mi">2</span><span class="p">)</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span> <span class="o">/</span> <span class="n">common_term</span><span class="p">)</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="p">(</span><span class="n">common_term</span><span class="p">)</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span> <span class="o">/</span> <span class="mi">12</span>
    <span class="n">t3</span> <span class="o">=</span> <span class="p">((</span><span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">+</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span> <span class="o">-</span> <span class="mi">2</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">loc</span><span class="p">)</span> <span class="o">/</span> <span class="mi">2</span><span class="p">)</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">t1</span> <span class="o">+</span> <span class="mf">0.5</span> <span class="o">*</span> <span class="p">(</span><span class="n">t2</span> <span class="o">+</span> <span class="n">t3</span><span class="p">)</span> <span class="o">/</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="mi">2</span><span class="p">)</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Uniform</span><span class="p">,</span> <span class="n">Pareto</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_uniform_pareto</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="n">support_uniform</span> <span class="o">=</span> <span class="n">p</span><span class="o">.</span><span class="n">high</span> <span class="o">-</span> <span class="n">p</span><span class="o">.</span><span class="n">low</span>
    <span class="n">t1</span> <span class="o">=</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">alpha</span> <span class="o">*</span> <span class="n">q</span><span class="o">.</span><span class="n">scale</span><span class="o">.</span><span class="n">pow</span><span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">alpha</span><span class="p">)</span> <span class="o">*</span> <span class="p">(</span><span class="n">support_uniform</span><span class="p">))</span><span class="o">.</span><span class="n">log</span><span class="p">()</span>
    <span class="n">t2</span> <span class="o">=</span> <span class="p">(</span><span class="n">_x_log_x</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">high</span><span class="p">)</span> <span class="o">-</span> <span class="n">_x_log_x</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">low</span><span class="p">)</span> <span class="o">-</span> <span class="n">support_uniform</span><span class="p">)</span> <span class="o">/</span> <span class="n">support_uniform</span>
    <span class="n">result</span> <span class="o">=</span> <span class="n">t2</span> <span class="o">*</span> <span class="p">(</span><span class="n">q</span><span class="o">.</span><span class="n">alpha</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span> <span class="o">-</span> <span class="n">t1</span>
    <span class="n">result</span><span class="p">[</span><span class="n">p</span><span class="o">.</span><span class="n">low</span> <span class="o">&lt;</span> <span class="n">q</span><span class="o">.</span><span class="n">support</span><span class="o">.</span><span class="n">lower_bound</span><span class="p">]</span> <span class="o">=</span> <span class="n">inf</span>
    <span class="k">return</span> <span class="n">result</span>


<span class="nd">@register_kl</span><span class="p">(</span><span class="n">Independent</span><span class="p">,</span> <span class="n">Independent</span><span class="p">)</span>
<span class="k">def</span> <span class="nf">_kl_independent_independent</span><span class="p">(</span><span class="n">p</span><span class="p">,</span> <span class="n">q</span><span class="p">):</span>
    <span class="k">if</span> <span class="n">p</span><span class="o">.</span><span class="n">reinterpreted_batch_ndims</span> <span class="o">!=</span> <span class="n">q</span><span class="o">.</span><span class="n">reinterpreted_batch_ndims</span><span class="p">:</span>
        <span class="k">raise</span> <span class="ne">NotImplementedError</span>
    <span class="n">result</span> <span class="o">=</span> <span class="n">kl_divergence</span><span class="p">(</span><span class="n">p</span><span class="o">.</span><span class="n">base_dist</span><span class="p">,</span> <span class="n">q</span><span class="o">.</span><span class="n">base_dist</span><span class="p">)</span>
    <span class="k">return</span> <span class="n">_sum_rightmost</span><span class="p">(</span><span class="n">result</span><span class="p">,</span> <span class="n">p</span><span class="o">.</span><span class="n">reinterpreted_batch_ndims</span><span class="p">)</span>
</pre></div>

             </article>
             
            </div>
            <footer>
  

  

    <hr>

  

  <div role="contentinfo">
    <p>
        &copy; Copyright 2019, Torch Contributors.

    </p>
  </div>
    
      <div>
        Built with <a href="http://sphinx-doc.org/">Sphinx</a> using a <a href="https://github.com/rtfd/sphinx_rtd_theme">theme</a> provided by <a href="https://readthedocs.org">Read the Docs</a>.
      </div>
     

</footer>

          </div>
        </div>

        <div class="pytorch-content-right" id="pytorch-content-right">
          <div class="pytorch-right-menu" id="pytorch-right-menu">
            <div class="pytorch-side-scroll" id="pytorch-side-scroll-right">
              
            </div>
          </div>
        </div>
      </section>
    </div>

  


  

     
       <script type="text/javascript" id="documentation_options" data-url_root="../../../" src="../../../_static/documentation_options.js"></script>
         <script src="../../../_static/jquery.js"></script>
         <script src="../../../_static/underscore.js"></script>
         <script src="../../../_static/doctools.js"></script>
         <script src="../../../_static/language_data.js"></script>
     

  

  <script type="text/javascript" src="../../../_static/js/vendor/popper.min.js"></script>
  <script type="text/javascript" src="../../../_static/js/vendor/bootstrap.min.js"></script>
  <script src="https://cdnjs.cloudflare.com/ajax/libs/list.js/1.5.0/list.min.js"></script>
  <script type="text/javascript" src="../../../_static/js/theme.js"></script>

  <script type="text/javascript">
      jQuery(function () {
          SphinxRtdTheme.Navigation.enable(true);
      });
  </script>
 
<script>
  (function(i,s,o,g,r,a,m){i['GoogleAnalyticsObject']=r;i[r]=i[r]||function(){
  (i[r].q=i[r].q||[]).push(arguments)},i[r].l=1*new Date();a=s.createElement(o),
  m=s.getElementsByTagName(o)[0];a.async=1;a.src=g;m.parentNode.insertBefore(a,m)
  })(window,document,'script','https://www.google-analytics.com/analytics.js','ga');

  ga('create', 'UA-90545585-1', 'auto');
  ga('send', 'pageview');

</script>

<script async src="https://www.googletagmanager.com/gtag/js?id=UA-117752657-2"></script>

<script>
  window.dataLayer = window.dataLayer || [];

  function gtag(){dataLayer.push(arguments);}

  gtag('js', new Date());
  gtag('config', 'UA-117752657-2');
</script>

<img height="1" width="1" style="border-style:none;" alt="" src="https://www.googleadservices.com/pagead/conversion/795629140/?label=txkmCPmdtosBENSssfsC&amp;guid=ON&amp;script=0"/>


  <!-- Begin Footer -->

  <div class="container-fluid docs-tutorials-resources" id="docs-tutorials-resources">
    <div class="container">
      <div class="row">
        <div class="col-md-4 text-center">
          <h2>Docs</h2>
          <p>Access comprehensive developer documentation for PyTorch</p>
          <a class="with-right-arrow" href="https://pytorch.org/docs/stable/index.html">View Docs</a>
        </div>

        <div class="col-md-4 text-center">
          <h2>Tutorials</h2>
          <p>Get in-depth tutorials for beginners and advanced developers</p>
          <a class="with-right-arrow" href="https://pytorch.org/tutorials">View Tutorials</a>
        </div>

        <div class="col-md-4 text-center">
          <h2>Resources</h2>
          <p>Find development resources and get your questions answered</p>
          <a class="with-right-arrow" href="https://pytorch.org/resources">View Resources</a>
        </div>
      </div>
    </div>
  </div>

  <footer class="site-footer">
    <div class="container footer-container">
      <div class="footer-logo-wrapper">
        <a href="https://pytorch.org/" class="footer-logo"></a>
      </div>

      <div class="footer-links-wrapper">
        <div class="footer-links-col">
          <ul>
            <li class="list-title"><a href="https://pytorch.org/">PyTorch</a></li>
            <li><a href="https://pytorch.org/get-started">Get Started</a></li>
            <li><a href="https://pytorch.org/features">Features</a></li>
            <li><a href="https://pytorch.org/ecosystem">Ecosystem</a></li>
            <li><a href="https://pytorch.org/blog/">Blog</a></li>
            <li><a href="https://github.com/pytorch/pytorch/blob/master/CONTRIBUTING.md">Contributing</a></li>
          </ul>
        </div>

        <div class="footer-links-col">
          <ul>
            <li class="list-title"><a href="https://pytorch.org/resources">Resources</a></li>
            <li><a href="https://pytorch.org/tutorials">Tutorials</a></li>
            <li><a href="https://pytorch.org/docs/stable/index.html">Docs</a></li>
            <li><a href="https://discuss.pytorch.org" target="_blank">Discuss</a></li>
            <li><a href="https://github.com/pytorch/pytorch/issues" target="_blank">Github Issues</a></li>
            <li><a href="https://pytorch.org/assets/brand-guidelines/PyTorch-Brand-Guidelines.pdf" target="_blank">Brand Guidelines</a></li>
          </ul>
        </div>

        <div class="footer-links-col follow-us-col">
          <ul>
            <li class="list-title">Stay Connected</li>
            <li>
              <div id="mc_embed_signup">
                <form
                  action="https://twitter.us14.list-manage.com/subscribe/post?u=75419c71fe0a935e53dfa4a3f&id=91d0dccd39"
                  method="post"
                  id="mc-embedded-subscribe-form"
                  name="mc-embedded-subscribe-form"
                  class="email-subscribe-form validate"
                  target="_blank"
                  novalidate>
                  <div id="mc_embed_signup_scroll" class="email-subscribe-form-fields-wrapper">
                    <div class="mc-field-group">
                      <label for="mce-EMAIL" style="display:none;">Email Address</label>
                      <input type="email" value="" name="EMAIL" class="required email" id="mce-EMAIL" placeholder="Email Address">
                    </div>

                    <div id="mce-responses" class="clear">
                      <div class="response" id="mce-error-response" style="display:none"></div>
                      <div class="response" id="mce-success-response" style="display:none"></div>
                    </div>    <!-- real people should not fill this in and expect good things - do not remove this or risk form bot signups-->

                    <div style="position: absolute; left: -5000px;" aria-hidden="true"><input type="text" name="b_75419c71fe0a935e53dfa4a3f_91d0dccd39" tabindex="-1" value=""></div>

                    <div class="clear">
                      <input type="submit" value="" name="subscribe" id="mc-embedded-subscribe" class="button email-subscribe-button">
                    </div>
                  </div>
                </form>
              </div>

            </li>
          </ul>

          <div class="footer-social-icons">
            <a href="https://www.facebook.com/pytorch" target="_blank" class="facebook"></a>
            <a href="https://twitter.com/pytorch" target="_blank" class="twitter"></a>
            <a href="https://www.youtube.com/pytorch" target="_blank" class="youtube"></a>
          </div>
        </div>
      </div>
    </div>
  </footer>

  <div class="cookie-banner-wrapper">
  <div class="container">
    <p class="gdpr-notice">To analyze traffic and optimize your experience, we serve cookies on this site. By clicking or navigating, you agree to allow our usage of cookies. As the current maintainers of this site, Facebook’s Cookies Policy applies. Learn more, including about available controls: <a href="https://www.facebook.com/policies/cookies/">Cookies Policy</a>.</p>
    <img class="close-button" src="../../../_static/images/pytorch-x.svg">
  </div>
</div>

  <!-- End Footer -->

  <!-- Begin Mobile Menu -->

  <div class="mobile-main-menu">
    <div class="container-fluid">
      <div class="container">
        <div class="mobile-main-menu-header-container">
          <a class="header-logo" href="https://pytorch.org/" aria-label="PyTorch"></a>
          <a class="main-menu-close-button" href="#" data-behavior="close-mobile-menu"></a>
        </div>
      </div>
    </div>

    <div class="mobile-main-menu-links-container">
      <div class="main-menu">
        <ul>
          <li>
            <a href="https://pytorch.org/get-started">Get Started</a>
          </li>

          <li>
            <a href="https://pytorch.org/features">Features</a>
          </li>

          <li>
            <a href="https://pytorch.org/ecosystem">Ecosystem</a>
          </li>

          <li>
            <a href="https://pytorch.org/mobile">Mobile</a>
          </li>

          <li>
            <a href="https://pytorch.org/hub">PyTorch Hub</a>
          </li>

          <li>
            <a href="https://pytorch.org/blog/">Blog</a>
          </li>

          <li>
            <a href="https://pytorch.org/tutorials">Tutorials</a>
          </li>

          <li class="active">
            <a href="https://pytorch.org/docs/stable/index.html">Docs</a>
          </li>

          <li>
            <a href="https://pytorch.org/resources">Resources</a>
          </li>

          <li>
            <a href="https://github.com/pytorch/pytorch">Github</a>
          </li>
        </ul>
      </div>
    </div>
  </div>

  <!-- End Mobile Menu -->

  <script type="text/javascript" src="../../../_static/js/vendor/anchor.min.js"></script>

  <script type="text/javascript">
    $(document).ready(function() {
      mobileMenu.bind();
      mobileTOC.bind();
      pytorchAnchors.bind();
      sideMenus.bind();
      scrollToAnchor.bind();
      highlightNavigation.bind();
      mainMenuDropdown.bind();
      filterTags.bind();

      // Remove any empty p tags that Sphinx adds
      $("[data-tags='null']").remove();

      // Add class to links that have code blocks, since we cannot create links in code blocks
      $("article.pytorch-article a span.pre").each(function(e) {
        $(this).closest("a").addClass("has-code");
      });
    })
  </script>
</body>
</html>